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ABSTRACT 

Context. The knowledge of the orbit or the ephemeris uncertainty of asteroid presents a particular interest for various purposes. These 
quantities are for instance useful for recovering asteroids, for identifying lost asteroids or for planning stellar occultation campaigns. 
They are also needed to estimate the close approach of Near-Earth asteroids, and subsequent risk of collision. Ephemeris accuracy 
can also be used for instrument calibration purposes or for scientific applications. 

Aims. Asteroid databases provide information about the uncertainty of the orbits allowing the measure of the quality of an orbit. The 
aims of this paper is to analyse these different uncertainty parameters and to estimate the impact of the different measurements on the 
uncertainty of orbits. 

Methods. We particularly deal with two main databases ASTORB and MPCORB providing uncertainty parameters for asteroid orbits. 
Statistical methods are used in order to estimate orbital uncertainty and compare with parameters from databases. Simulations are also 
generated to deal with specific measurements such as future Gaia or present radar measurements. 

Results. Relations between the uncertainty parameter and the characteristics of the asteroid (orbital arc, absolute magnitude, ...) are 
highlighted. Moreover, a review of the different measuments are compiled and the impact of these measures on the accuracy of the 
orbit is also estimated. 

Key words. Minor planets, asteroids: general - Ephemerides - Astrometry 



1. Introduction 

Asteroids, since the time of their first discoveries in the nine- 
teenth century, have raised the need to compute and to predict 
their position for ensuring subsequent detections and observa- 
tions. While methods for computing orbits or predicting posi- 
tions of planets and comets had been developed at that time, 
the case of asteroids (or minor planets) raised new challenges; 
because no specific assumption can be made on the orbital ele- 
ments, that are all completely unknown. This was successfully 



solved by GauB (Gauss 1809, 1864) and his method of orbit 



computation from the knowledge of three topocentric positions. 
That method, presented as "a judicious balancing of geometri- 
cal and dynamical concepts" by Escobal (1965 ), proved to be 
remarkably efficient providing — shortly after the discovery of 
Ceres by Piazzi, in 1801, and observations over only a 40 days 
covering a 3 degree arc — a prediction for the next apparition of 
Ceres twelve month later, to better than 0.1 degree. Without be- 
ing exhaustive, one can also mention here, for completion, the 
method of Laplace (e.g. Poincare 1906 ), based on the knowledge 
of a position and its derivatives, which shows several similarities 
on the fundamental scheme (Tisserand & Perchot 1899; Celletti 
& Pinzari|2005]). These methods ([Herget|1948||Dubyago|196r 



Escobal|1965HRoy|2005||Milani & Gro nchi 20 10J and reference 
therein) paved the way for many more discovery of minor plan- 
ets, or asteroids, and continuous observations of these bodies. 
The number of known asteroids is still increasing, presently at 
the end 2012 there are slightly less than 600,000 catalogued bod- 
ies in the astorb database ( |Bowell|2012| ); and it is interesting to 
remind the remark extracted from the preface of E. Dubois, the 



French translator of Gauss' work (Gauss 1864]): 
"Or il est bien probable que la zone situee entre Mars et Jupiter 
n'est pas encore suffisamment exploree et que le chiffre de 79 
auquel on est arrive, sera encore augmente. Qui sait ce que 
reserve l'avenir !!... Bientot alors les astronomes officiels n'y 
pourront plus suffire, si des calculateurs devoues a l'astronomie 
et a ses progres ne leur viennent aussi en aide de ce coteQ'. 
Preliminary orbit computations are thus needed to provide some 
initial knowledge of the celestial object's orbit. They differ in 
nature, as an inverse problem, from the direct problem of subse- 
quent position prediction and orbit propagation that can be more 
complete in terms of correction to the observations and force 
models. 

Nowadays, ephemerides are commonly needed for various 
practical or scientific use: to prepare an observation program or 
space mission rendez-vous, to be able to cross-match a source 
or to identify a known asteroid observed in a CCD frame, to 
link two observed arcs as being of the same object, for predict- 
ing stellar occultations, for developing planetary ephemerides, 
for testing new dynamical models and physics, for contribut- 
ing to local tests of the general relativity, for deriving asteroid 
masses from close encounters, or non-gravitationnal parameters, 
etc. During previous centuries these objects were considered as 
planets or minor planets, and analytical or semi-numerical theo- 
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1 "It is likely that the zone between Mars and Jupiter is not yet suf- 
ficiently explored and that the number of 79 that has been reached will 
still be increased. Who knows what the future will hold for us!! Soon 
official astronomers will not suffice, if some calculators devoted to as- 
tronomy and its progress do not come to their rescue also on that side." It 
is clear that two centuries later modern electronic computing machines 
came to the rescue and are now impossible to circumvent. 
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ries could be applied to each body before the advent of modern 
powerful computers (e.g. Leveau|1880; Brouwer 1937 Morando 



1965). Ephemerides are now computed directly by numerical in- 
tegration of the equation of motions for the orbit propagation us- 
ing different integrators. Such numerical integration requires the 
use of a standard dynamical model and associated parameters for 
all solutions to be consistent, and the knowledge of initial condi- 
tions to the Cauchy problem as given by orbital parameters at a 
reference epoch. The latter are obtained, as for any solar system 
object, by a fit to the observations and are updated regularly. The 
terminology 'ephemerides' is broadly used for any quantity re- 
lated to the computation of such orbit propagation of a celestial 
object, giving its position and motion (spherical cartesian, appar- 
ent, astrometric, as seen from various centre, in various reference 
planes, at various dates, on different time scales, etc.). 

In addition to these calculated values from orbit propaga- 
tion and subsequent transformations corresponding generally to 
a nominal solution (e.g. least-squares with a standard dynamical 
model), knowledge of the ephemeris uncertainty, precision and 
accuracy is often mandatory as is the case for any physical quan- 
tities. It is also often useful to know the degree of confidence one 
can have in a predicted quantity. Further, the ephemerides uncer- 
tainty is needed in the various fields presented above and in par- 
ticular for planning observations using instruments with small 
field of view, for instrument calibration purposes, or when one 
needs to quantify the probability of an impact of a Near-Earth 
Object (NEO) with a terrestrial planet, or for efficient planning 
of stellar occultation campaigns (in particular if it involves large 
number of participants and/or large telescopes), etc. 

The precision of an ephemeris can be incorrectly reduced to 
the precision of the observations it is based on. Of course, all 
things being equal, the higher the precision of the observation 
the better the theory and ephemerides; but this will not yield an 
indication of how fast the precision for any predicted or extrap- 
olated quantity is being degraded. Indeed, the precision of the 
ephemerides is a quantity varying with time; when extrapolat- 
ing the motion to dates far into the future or in the past (several 
centuries), the precision will globally be poorer; up to the point 
where chaotic orbits on time span of million years make such 
ephemeris prediction unrealistic. An ephemeris uncertainty can 
be decreased on a short time span — for instance if required by 
a flyby or space mission to a solar system body — with last min- 
utes observations. However, if no more observations are added 
in the fitting process, the ephemerides precision will inevitably 
degrade in time. The accuracy of an ephemeris as the precision 
of the theory for the dynamical model and its representation can 
be internal or external. Internal precision refers to the numerical 
model used (terms of the development in an analytical theory or 
machine precision in case of numerical integration of the ODE 
for the equations of motions) or additionally to the representa- 
tion used to compress such ephemerides (Chebychev polynomi- 
als, mixed functions, etc.). External precision refers to the ad- 
justment to the observations and hence to the stochastic and sys- 
tematic errors involved in the measurements and data reduction, 
in the validity of the models used to fit and transform the data, in 
the uncertainty in the considered parameters, etc. We will deal in 
the following with the external errors, which accounts for most 
of the uncertainty in present ephemerides of asteroids and other 
dwarf planets, and small bodies of the Solar System. 

All ephemerides require observations and measurents related 
to the dynamical model of the equations of motion. Current 
ground-based surveys (LINEAR, Catalina, Spacewatch, C9) 
which where basically designed to detect 90% of the largest 
NEOs provide most of such data; they are completed by some 



scientific programs such as radar observations from Arecibo 
and Goldstone for NEOs the CFHTLS Ec liptic Survey for the 
Trans -Neptunian Objects (Jones et al. 2006 ) or the Deep Ecliptic 
Survey (DES) for Kuiper Belt Objects and Centaurs ( |Elliot et aT] 
2005[ ). In the very near future surveys such as Pan-starrs, Gaia, 
LSST, will also provide large amount of astrometric positions. 

In this paper, we investigate the orbital uncertainty of aster- 
oids in the numerical and statistical way. First, we briefly present 
the two main databases of orbital parameters of asteroids. In par- 
ticular, we present and compare uncertainty parameters provided 
by these databases (Sect. [2]l. 

In Sect. [3] we highlight relations between orbital uncertainty 
and other asteroid parameters such as dynamical classes of as- 
teroids, magnitude. We specifically study the relation between 
orbital arc and uncertainty, in particular for short orbital arcs. 
Eventually, in Sect. [4] we present statistical information about 
astrometric measurements of asteroids and we quantify the im- 
pact of astrometric measurements in the radar and Gaia space 
mission context. 



2. Ephemeris uncertainty parameters 

The number of discovered asteroids has exceeded 590,000 at 
the beginning of October 2012, and the discovery rate is still 
of about 1800 new asteroids per montrj^] Currently, four main 
centres provide asteroid orbital databases (Minor Planet Center, 
Lowell Observatory, Jet Propulsion Laboratory, Pisa University). 
In this paper, we mainly deal with two of these databases: as- 
torb from Lowell Observatory ( |Bowell|2012| i and mpcorb from 
IMinor Planet Centerl (|2012bb. A total of 590,095 asteroids are 



compiled in astorb and 590,073 in mpcorb as from October 5th 
2012. The two databases provide similar asteroid parameters, in 
particular: 

- name/number: name or preliminary designation, asteroid 
number; 

- osculating elements: semi-major axis a, eccentricity e, incli- 
nation /, mean anomaly M, argument of perihelion to, longi- 
tude of ascending node Q at reference epoch; 

- magnitude: absolute magnitude H, slope parameter G; 

- observations: number of astrometric observations, arc 
length or year of first and last observation. 

The two databases also provide parameters about pre- 
dictability of the ephemerides. Only one for mpcorb which is 
the U parameter and five parameters for astorb which are CEU, 
PEU, next PEU and the two greatest PEU. 

The uncertainty parameter U is an integer value between 
and 9, where indicates a very small uncertainty and 9 an 
extremely large uncertainty. Detailed information about the 
computation is given in Mi nor Planet Center] ( |2012d| l. Briefly, 
the U parameter is computed thanks to another parameter, 
RUNOFF, which depends on the uncertainty in the time of 
passage at perihelion, the orbital period and its uncertainty. 
RUNOFF expresses the uncertainty in longitude in seconds of 
arc per decade. The U parameter is derived by a logarithmic 
relation of RUNOFF. The quality of the orbit can be quickly 
determined with the uncertainty parameter U. 

The astorb database provides five parameters related to the 
ephemeris uncertainty. In our study, we specifically dealt with 
two of them, current ephemeris uncertainty (CEU) and the rate 



Average of the first half of year 2012. 



2 



J. Desmars et al.: Statistical and Numerical Study of Asteroid Orbital Uncertainty 



of change of CEU. The CEU matches to the sky-plane uncer- 
tainty cr^ at a daten provided by astorb. A brief description of 
the uncertainty-analysis technique is presented in|Yeomans et al. 



1987) and all the details can be found in IMuinonen & Bowell 



1993). The orbit determination provides the covariance matrix 



of the orbital elements A. Linear transformations then give the 
covariance matrix X in spherical coordinates (in right ascension 
and declination): 



2 = ¥A*F 



(1) 



where *P is the matrix of partial derivatives between spherical 
coordinates and orbital elements. Finally, by propagating this co- 
variance matrix at a given date (date of CEU), the sky-plane un- 
certainty can be determined as the trace of the matrix ( Mui nonen] 
& Bowell|1993[l. The CEU is determined in two body linear ap- 
proximation ( Bowell||2012 1. The rate of change of CEU (noted 
CEU rate) is in arcsec/day. According to Muin onen & BoweII| 
( 1993 ), in linear approximation, uncertainties in spherical coor- 
dinates increase linearly with time in the two body approxima- 
tion. 



2.1. Comparison of the ephemeris uncertainty parameters 

The U parameter can be compared to the CEU and its rate of 
change. Figure [T] shows a correlation between these parameters. 
Qualitatively and as expected, asteroids with a low U parameter 
have also a low CEU and a CEU rate. 
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Fig. 1. CEU, CEU rate, U parameter for asteroids in common to 
the astorb and mpcorb databases. 



The measure of the uncertainty could be provided by only 
one of these parameters. Nevertheless, those two parameters are 
not perfect and some problems can appear. 

Indeed, the main difficulty with the uncertainty parameters 
is that they are sometimes not indicated. For mpcorb, 105,399 
asteroids (that is 17.9 % of the total in the database) have no U 
parameter at all and 9431 (1.6 %) have only a qualitative let- 
ter for the U paramete^] In the astorb database, the CEU has 

3 The date of CEU is usually from to 40 days before the epoch of 
osculation depending on the update of the database. In our case, the date 
of CEU and the epoch of osculation are the same September 30th 2012. 
For previous updates, the difference can reach about 45 days 

4 In the mpcorb database, the U parameter can be indicated as a let- 
ter. If U is indicated as 'E', it means that the orbital eccentricity was 



not been computed (and appears as 0) for 2585 asteroids (0.44 
%). Moreover, CEU is determined using the approximate two- 
body problem. For some critical cases (Earth-approaching aster- 
oids), the uncertainties may have been misestimated by a fac- 
tor of several (Bowell 2012). As an illustration, we have com- 
puted the orbital uncertainty - with methods described in the 
next section - for the asteroid Apophis which is a well known 
Earth approaching PHA. The difference between the simplified 
two-body approximation and the full N-body perturbations is 
clear (see Fig. EEl, The ratio between the two values of uncer- 
tainty is close to 1 until 2029 and reaches approximately 10 5 
after the 2029-close approach. The 2-body approach does not 
involve the stretching of the ortbital uncertainty and thus yield 
to over-optimistic results. 
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Fig. 2. Orbital uncertainty in distance of asteroid Apophis using 
two-body approximation and full N-body perturbations. 



2.2. Study of the ephemeris uncertainty parameters 

Another problem is that some asteroids have CEU and U param- 
eter in total contradiction. For example, they can have a good 
CEU (less than 1 arcsec) and a bad U parameter (U=9) (and 
conversly). To clarify this situation, we have compared the CEU 
value to the standard deviation <x obtained by orbital clones for 
eight selected representative asteroids (2 with a bad U and bad 
CEU, 2 with bad U and good CEU, 2 with good U and bad CEU 
and 2 with good U and good CEU). 

In this context, we have computed the standard deviation 
provided by clones with two different methods using non lin- 
ear extrapolation. The main principle of these methods is to per- 
form a Monte Carlo propagation of the orbit, i.e. to compute 
clones of initial conditions around the nominal orbit, providing 
as many orbits as possible. The first method (Monte Carlo using 
Covariance Matrix, MCCM) consists in adding a random noise 
to the set of nominal initial conditions using the standard devia- 
tion and correlation of these parameters thanks to the covariance 
matrix. The second one (Bootstrap Resampling, BR) uses the 
bootstrap resampling method directly on the observations, and 



assumed. For one-opposition orbits, U can also be 'D' if a double (or 
multiple) designation is involved or 'F' if an assumed double (or multi- 
ple) designation is involved (Minor Planet Center 201 2d} . 
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consists in determining new set of initial conditions by fitting to 
a bootstrap set of observations. These two methods, one para- 
metric and one non parametric, have been used previously in 
|Desmars et al.| ( p009] l for the study of the precision of planetary 
satellites ephemerides. 

For this test, we have computed 1000 clones for orbital initial 
conditions of each representative asteroid. Then we have calcu- 
lated the standard deviation <x of the angular separation s in the 
plane-of-sky at the date of CEU, which represents the angular 
deviation of an orbit to the nominal one, defined as: 



si = V((Q'; - »o) cos Si) 2 + {Si - S ) 2 



1 



(2) 



where i denotes the number of the clone, the subscript 
refers to the nominal orbit, s is the mean, and N is the total 
number of clones computed. 

Table[T]gives the comparison between the U parameter, CEU 
and standard deviation provided by these orbital clones for the 
two methods. 



Table 1. Comparison of the U parameter, the CEU and stan- 
dard deviation <x obtained from orbital clones around the nomi- 
nal orbits, with two different methods (MCCM and BR, see text). 
Units are in arcseconds for CEU and cr, U is unitless. 





U 


CEU 


°" MCCM 


cr B R 




(mpcorb) 


(astorb) 


(our work) 


(our work) 


2002GM5 


9 


180,000 


161,461 


122,517 


2006LA 


9 


340,000 


211,970 


235,317 


2000RH60 


5 


0.17 


0.216 


0.275 


2010JK1 


6 


0.68 


0.451 


0.367 


2007WD5 





4300 


67,388 


14,844 


2003TO9 





240 


161 


128 


4321 Zero 





0.058 


0.077 


0.086 


31824 Elatus 





0.36 


0.548 


0.331 



For these representative asteroids, the CEU is often close to 
the standard deviation computed here, whereas the U param- 
eter is misestimated for at least four representative asteroids. 
Moreover, we stress that the U parameter provides a number 
which is not related to a physical value (a distance or an angle). 

In light of the previous tests for several particular cases, and 
of Fig{T]for the general purpose, the CEU seems to be a good and 
practical parameter to estimate the accuracy and predictability 
of an asteroid orbit, as CEU is quickly computable, precise and 
providing a physical value (an angle). 



3. Relations between CEU and physical or orbital 
parameters 

As the CEU seems to be a useful parameter to measure the or- 
bital uncertainty, we will hereafter use this parameter and high- 
light several relations between CEU and orbital parameters. 



3.1. Relation between absolute magnitude, orbital arc and 
CEU 

Generally, asteroids can be classified in dynamical groups. We 
propose, for sake of simplicity, to gather the asteroids into the 
following three main groups defined as^] 

- A-NEA: asteroid with perihelion q < 1 .3 au. It represents 
asteroids in inner zone of solar system. This group includes 
the Aten, Apollo, Amor, and Inner Earth Orbit (IEO). There 
are 9144 objects in A-NEA group in astorb. In this cate- 
gory, Potentially Hazardous Asteroids (PHA) are asteroids 
presenting a risk of devastating collision (i.e. with H < 22 
and a Minimum Orbital Intersection Distance, MOID< 0.05 
au). 

- B-MBA: asteroid with perihelion q > 1 .3 au and semi-major 
axis a < 5.5 au, representing asteroids in intermediate zone. 
Thus Trojans are considered as B-MBA in the following. 
There are 579,208 B-MBAs in astorb. 

- C-TNO: asteroid with perihelion q > 1.3 au and semi -major 
axis a > 5.5 au, representing asteroids in outer zone of so- 
lar system. Thus Centaurs, asteroids between Jupiter and 
Neptune, are included in C-TNO. There are 1743 C-TNOs 
in astorb. 

Figure [3]represents the absolute magnitude H, CEU and the 
length of orbital arc given in astorb for all the asteroids. Each 
group defined previously appears to be clearly distinguishable: 
C-TNOs in the left part (cross), B-MBAs in middle part (small 
square) and A-NEAs in right part (bullet). Large C-TNOs (larger 
than »45 km corresponding to absolute magnitude 8, depending 
on albedo) with short arcs and small A-NEAs (approximately 
smaller than 5 km whatever the albedo of the asteroid larger 
than 0.05, using the relation between magnitude and size (Minor 
Planet Center|2012c i) with short arcs also appear. 
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Fig. 3. H magnitude, CEU and orbital arc for asteroid from as- 
torb. Cross is for C-TNO, small square for B-MBA and bullet 
for A-NEA. When H or CEU are unknown, they are represented 
on the plot with the values H — 35 or CEU= 10 9 , in contrast to 
the convention adopted in astorb database where H = 14. 



One also notes, that for all the objects, a correlation between 
the value of CEU, the length of orbital arc and the number of 
observations can be identified. In Figj4] this correlation appears 



Number of objects are given at the date of October 5th 2012. 
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clearly. We can identify three groups related to the length of the 3.2. Relation between orbital arc and CEU 



arc < 10 days: these asteroids are barely observed (less 
than 10 observations), corresponding mostly to the discov- 
ery period, and have a high CEU (between 1 and 10 7 arcsec). 
Moreover, the CEU does not decrease much on such graph 
with increasing orbital arc (see Sect 3.3 for a discussion of 
this statement); 

10 days < arc < 250 days: these asteroids are not much 
observed (less than 100 observations) and have a high CEU 
(between 1 and 10 6 arcsec); 

250 days < arc < 8000 days: these asteroids are also not 
much observed but have a better CEU (between 0.1 and 100 
arcsec); 

8000 days < arc: these asteroids have all a large number of 
observations (more than 100 observations) and have a good 
CEU (between 0.01 and 10 arcsec). 
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Fig. 4. Relation between the length of the orbital arc, CEU and 
the number of observations for all asteroids. 



The gap in orbital arc apparent around 250 days means that 
when an asteroid is discovered, it is rarely observed for more 
than 250 days. Beyond this period, the asteroid is in unfavorable 
configuration to be observed because of its small solar elonga- 
tion (conjunction) and can be observed again only after its pas- 
sage behind the Sun. Two peaks can also be observed in Fig(4] 
at orbital arc longer than 10,000 days. One for arcs of 19,000 
days and another, more spread out, for approximately 14,000 
days. These peaks can be explained by an important number 
of objects with such an orbital arc. Actually, many asteroids 
(~ 1500) have been discovered in 1960 with the 1.22m tele- 
scope of Mt.Palomar, similarly in 1973. Because these asteroids 
have been observed recently (in the late 2000s), their orbital arc 
correspond approximatively to 14,000 days (i.e. the period 1973- 
2012) and to 19,000 days (i.e. the period 1960-2012). These ob- 
jects have now various CEU corresponding to the peaks. 

With these considerations, asteroids can be classified in 12 
groups considering their orbital arc, their dynamical classes (A- 
NEA, B-MBA, C-TNO) and their CEU (given the correlation 
between CEU and length of arc). Table [2] gives the number and 
the median^] CEU for each group. Asteroids with undetermined 
CEU (about 2585 asteroids with CEU=0 in astorb) are not taken 
into account. 



Figure [4] shows that the decrease of the CEU is correlated to the 
increase of the orbital arc. Moreover the increasing rate of the 
CEU appears different with varying orbital arc. We have com- 
puted an empirical relation for this improvement rate for each 
of the groups defined in the previous section. For each group, a 
linear regression can be computed between CEU and the length 
of orbital arc (as shown in Pig j5j and we have the relation: 

log(C EU) = a logiarc) + b (3) 

where a and b are given in Tableland Fig|5]for each groups. 
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Fig. 5. Relation between the length of the orbital arc and CEU 
for all asteroids, and slope for the 4 classes identified (red for 
arc< lOd, green for lOd < arc < 250d, blue for 250d < arc < 
8000d and orange for arc > 8000d). 



If the orbital arc is smaller than 10 daysf 7 ] the CEU is 
very large and does not improve much when the length of arc 
becomes larger (see also Sect |3.3 I. However, in such a case, 
observations remain useful to avoid loosing the object. Between 
10 and 250 day arcs, the CEU is clearly improved when the 
length of arc increases. For asteroids with an orbital arc between 
250 and 8000 days, the CEU is smaller than 100 arcsec and is 
still much improved when the length of arc becomes greater. If 
the orbital arc is greater than 8000 days, then the CEU is not 
much improved and reaches its typical minimum value of about 
0.1-0.2 arcsec. 

As noted previously, in the short arc group (orbital arc less 
than 10 days) the CEU does not decrease quickly while the or- 
bital arc increases. It is unclear however if this is associated to 
the population on average or on each individual arc. In the latter 
case it could be a surprising result; so we highlight the trend on 
specific study case in the following section. 

3.3. Evolution of CEU for short orbital arc 

In this section, we investigate the evolution of the CEU as a func- 
tion of the orbital arc, for a single object. We try to highlight 



6 The median is to be prefered to the mean that can be - and generally 
is - dominated by few large values. 



Orbital arcs indicated as in the database are considered as 0. 1 in 
this section to allow the use of logarithmic scale. 
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Table 2. Different classes of ephemerides precision as given by their median current ephemeris uncertainty (CEU) and length of 
arc. Numbers of each category are in subscript. Asteroids with undetermined CEU are not taken into account here. 





B-MBA[57 6j6 47] 


A-NEA [9123] 


C-TNO [1711] 


arc< 10 [5 o,908] 
10 <arc< 250[8i,54O] 
250 <arc< 8000 [4 i 8 ,33i] 
8000 < arc [36 ,73i] 


380,000(48,326] 
1 100[77,579] 
0.11 [414,410] 
0.045(36,332] 


6300 [2 ,33] 

530[3587] 
0.16 B 95i] 
0.044 [35 2] 


420,000 [34 9] 
370,000(374] 

1 -0(970] 
0.13(47] 



Table 3. Values of the linear regression \og{CEU) - a log(arc)- 
b, for each groups defined in Tabl^2] 



arc (days) 


a 


b 


number 


0.1-10 


-1.448 


+5.804 


50,908 


10-250 


-2.460 


+6.774 


81,540 


250-8000 


-2.036 


+6.528 


418,331 


> 8000 


-0.254 


-0.270 


36,731 



the time spent after the discovery to really improve the quality 
of the orbit. A previous work in |Virtanen & Muinonen ([2006 ) 
studied the time evolution of the orbital elements uncertainty of 
the asteroid 2004 AS 1 together with the evolution of the impact 
risk assessment. For this purpose, the authors used the volume 
of variation (VoV) technique which is a nonlinear Monte Carlo 
technique (Muinone n et al.|2006[ ). 

In this context, we computed the orbital uncertainty of an 
asteroid by considering a variable (and short) length of orbital 
arc. For our work, six test case objects were considered: (31824) 
Elatus and (2060) Chiron for Centaur asteroids; (1) Ceres and 
(8) Flora for MBAs; (1866) Sisyphus and (4179) Toutatis (PHA) 
for NEOs. First, we performed a preliminary orbit determination 
using the Gauss method. The Gauss method converged two days 
after discovery for all cases except for Chiron and Ceres which 
converged one day after discovery. To assess the uncertainty of 
the orbit, we used the differential correction method. All of these 
computations were performed using OrbFit4.2 packag^] Then, 
for each arc, we computed 1001 clones of the nominal orbit from 
a multivariate Gaussian distribution given by the covariance ma- 
trix (obtained from the procedure described above). Each clone 
orbit is further directly propagated at the date of computation of 
CEU with Monte Carlo technique. The Lie integrator (Ba nceliri] 
et al. 2012c I was used to perform this propagation. Thus, we 
were able to compute the standard deviation cr using (Eq. [2J 
that we assumed to match with the definition of the CEU. If the 
Gauss method and/or orbital improvement failed or if the cr value 
obtained is large (greater that 10 7 ), then the CEU is considered 
as indeterminate and is set to 10 7 . Obviously, some asteroids 
considered in this study have an indeterminate CEU at the epoch 
of their discovery (Fig.|6]l. For Chiron and Ceres, a preliminary 
solution was obtained but the uncertainty on the orbit was too 
large. For the NEOs, the Gauss method converged 2 days after 
the discovery and the differential correction succeeded. For the 
MBAs, it is defined after 3 days whereas for the Centaurs, the 
uncertainty is not defined up to ten days. Then, the orbital uncer- 
tainty drops almost linearly for all the objects: for the Centaur 
asteroids, it goes below 100-1000 arcsec for an arc length 800- 
1000 days. This is due to the lack of observations because they 
are distant objects and they have a slow apparent motion. For 



the other asteroids, the orbital uncertainty decreases linearly and 
goes below 100 arcsec. 

One can see that its variation for the asteroid Toutatis has a 
stranger behaviour than the others because the sharp drop of the 
uncertainty value occurs around the 14th day. 



(31824) Elatus [Centaur] 
(2060) Chiron [Centaur] 
(8) Flora [MBA] 
(1) Ceres [MBA] 
(1866) Sisyphus [NEO] 
(4179) Toutatis [NEO, PHA] 



~ 100000 




http : / /adams . dm . unipi . i t/orbf it/ 



arc (days) 

Fig. 6. Variation of the orbital uncertainty as a function of the 
arc length for six single objects. The axis are expressed in a log- 
arithmic scale. 

This test shows that, except for the distant objects (Centaurs), 
observations performed after their discovery (until ~ 6 days) will 
mostly avoid the loss of the object than improve its orbit, as the 
rate of change is not that significant. However, continuous ob- 
servations performed after ~ 6 days will allow a significant orbit 
improvement as the rate is much steeper. 



4. Astrometric measurements 

In this section, we first described the different types of measure- 
ment of asteroid position, then we estimate the effect of specific 
measurements (optical observations, radar or future Gaia obser- 
vations) in the improvement of the orbit. 

4.1. Presentation 

The position of an asteroid can be determined by many types of 
measurements. The Minor Planet Center provides all the mea- 
surements related to asteroids for numbered and unnumbered 
asteroids (Minor Planet Center 2012a I. As of September 30th 
2012, 94,082,128 observations are available on the database. 
Table |4] gives the statistics about the different kind of measure- 
ments provided by the MPC database in particular their percent- 
age and the time span. 
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Historically, the first observations of asteroids were made 
with micrometer during the 19 th century. In the latter part of that 
century, the photographic plates allowed a better measurement of 
asteroid position. In the mid-1980s, the CCD has revolutionized 
astrometric measurement thanks to a better sensitivity, stability, 
dynamic range and ease of extracting astrometric measurements. 
Consequently, most of the measurement in the MPC database 
come from CCD. Space observations such as by Hubble Space 
Telescope, Spitzer and particularly WISE also provide an impor- 
tant number of measurements. 

In addition to the different types of measurements mentioned 
above, some observations can appear as secondary or marginal 
for different reasons; even if being precise they can have a lower 
impact on the orbital models. For example, the space mission 
Hipparcos also provides geocentric data ( |Hestroffer et al.|1 998 ). 
Radar measurements are not very common but give the most 
accurate measurements for asteroid. They also can provide 
physical characteristics such as shape. Astrometric positions 
can also be deduced from stellar occultations. 

The accuracy of each measurement can be obtained thanks 



to the AstDyS website ( |AstDyS 2012 1. For each numbered aster- 
oid, this website provides an observation file indicating the resid- 
ual (O-C) of each measurement. All the observations of all num- 
bered asteroid^] have been compiled and the accuracy of each 
measurement has been estimated by computing a weightecp"] 
root mean square of all residuals of the measurement. 

Table [5] provides the number, the percentage of accepted 
measurement and the estimated accuracy for each kind of ob- 
servations. After radar measurement, observations derived from 
occultations and Hipparcos geocentric observations appear as 
the more accurate one (about 0.08-0.15 arcsec). Recent mea- 
surements such as CCD observations or meridian circle obser- 
vations are rather accurate (between 0.2-0.4 arcsec). Older mea- 
surements such as photographic or micrometric measures are 
less accurate (1.0-1.8 arcsec). Surprisingly, observations from 
spacecraft (HST, Spitzer or WISE) are not any more accurate 
than groundbased data. Nevertheless, Wide-field Infrared Survey 
Explorer (WISE) was designed for asteroid detection and actu- 
ally, the mission has detected 157,000 asteroids, including more 
than 500 NEO ( jMainzer et al.|2011) . 

4.2. Impact of astrometric measurements on orbit 
uncertainty - The Apophis example 

In order to quantify the impact of a type of measurement on 
the orbital uncertainty of an asteroid, we specifically dealt with 
asteroid (99942) Apophis. This particular object belongs to the 
PHAs family and is known to be the most threatening object of 
the last decade since it is the only asteroid who reached level 4 
of Torino scale for a potential impact with the Earth on April 
2029 (Sansaturio & Arratia 2008). Since, new observations 
(optical and radar) ruled out every possibility of impact in 2029 
and this threat turned out to be a close encounter within 5.64R ffi 
with Earth. But this close encounter is so deep that the asteroid 
will move from the Aten family to the Apollo family orbit. 
Besides, its orbit will become chaotic and new possibilities 
of collision with the Earth after 2029 will appear. The most 
consensual collision date is in 2036 for which, at the date of 




2012 2014 2016 2018 2020 2022 2024 2026 2028 
Time [years] 

Fig. 7. Position uncertainty on the geocentric distance of 
Apophis after fitting to optical observations (solid line) and op- 
tical and radar data observations (dashed line). 



the last observations available in 2012, the risk was estimated 
with a probability of ~ 10~ 6 . This date is also famous because 
of the size of keyhole, region in the B-plane (Vals ecchi et al.| 
2003 ) of 2029 where the asteroid has to pass in order to collide 
the Earth in the future. The 2036-keyhole size is estimated at 
~ 600m. In the B-plane, we cannot only represent the state of 
the asteroid, located by two geocentric coordinates (£, £), at 
the date of close encounter, but also its uncertainty (cr^, cr^) 
and the relative position of the keyhole with respect to the 
nominal solution. More often than not, the uncertainty region 
is a 3cr-ellipse centered on the nominal solution and its size is 
directly linked to the observations used. Therefore, the position 
of keyhole compared with the size of the ellipse uncertainty is 
also important to quantify the risk of future collision (Chesley 
[20061 IKochetova et al.|2009l|Bancelin et al.|2012a) . 



9 As from October 5th 2012, there are 337,008 numbered asteroids 
in AstDyS database. 

10 AstDys provides a weight for each observation according to their 
accuracy. 



4.2.1. Past observations 

In this section, we first deal with the impact of the observations 
of Apophis (from 2004-2008), as available at MPC-UAI, on the 
uncertainty region. We propagated the equations of motion (as 
well as the variational equations) of Apophis until 2029, consid- 
ering gravitational perturbations of all the planets, including also 
the gravitational perturbations of the Moon, Ceres, Pallas and 
Vesta. Relativistic accelerations were also taken into account. 
The uncertainty of the nominal orbit was propagated using a lin- 
ear propagation of the initial covariance matrix in order to assess 
the evolution of the orbit uncertainty. We performed this test us- 
ing two different sets of observations. The first one uses only op- 
tical data from 2004-2008 and the second one uses optical and 
radar data (5 radar measurements - ranging and Doppler mea- 
surements - were performed at Arecibo in 2005 and 2006). This 
test is done in order to highlight the impact of radar data on the 
current orbit uncertainty. 

First, we represent in Fig. [7] the current position uncertainty of 
Apophis until the date of close encounter, when considering op- 
tical data only or optical and radar data. One can see that the 
uncertainty is improved by a factor of more than 2 when consid- 
ering radar data. This is not a surprise in as much as, radar data 
are very accurate measurements. 
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Table 4. Statistics of different kind of measurement provided by MPC database for numbered and unnumbered asteroids. 



code 


type 


number 


percentage 


timespan 


C 


CCD 


88,546,921 


94.12% 


1986-2012 


S/s 


Space observation 


4,006,572 


4.26% 


1994-2011 




HST 


3544 


0.00% 


1994-2010 




Spitzer 


114 


0.00% 


2004-2004 




WISE 


4,002,914 


4.25% 


2010-2011 


A 


Observations from B 1950.0 converted to J2000.0 


647,690 


0.69% 


1802-1999 


c 


Corrected without republication CCD observation 


462,065 


0.49% 


1991-2007 


P 


Photographic 


352,449 


0.37% 


1898-2012 


T 


Meridian or transit circle 


26,968 


0.03% 


1984-2005 


A/X 


Discovery observation 


1 C HA 1 

16, /41 


n now 
0.02% 


i o a c on 1 n 
1845-zUlO 


M 


Micrometer 


12,081 


0.01% 


1845-1954 


H 


Hipparcos geocentric observation 


5,494 


0.01% 


1989-1993 


R 


Radar observation 


2,901 


0.00% 


1968-2012 


E 


Occultations derived observation 


1,737 


0.00% 


1961-2012 


V 


"Roving observer" observation 


388 


0.00% 


2000-2012 


n 


Mini-normal place derived from averaging observations from video frames 


105 


0.00% 


2009-2012 


e 


Encoder 


16 


0.00% 


1993-1995 



Table 5. Accuracy of measurements for numbered asteroids from AstDyS2 database. 



type 


name 


number of 
measurement 


percentage of 
accepted measurement 


accuracy 


C 


CCD 


79,569,190 


99.49% 


0.388 


arcsec 


S 


Wise 


1,526,466 


99.86% 


0.583 


arcsec 


s 


Hubble Space Telesc. 


867 


96.54% 


0.585 


arcsec 


s 


Spitzer 


48 


33.33% 


1.673 


arcsec 


A 


B 1950 to J2000 


632,428 


82.17% 


1.170 


arcsec 


c 


Corrected CCD obs. 


423,792 


99.70% 


0.507 


arcsec 


P 


Photographic 


346,947 


93.38% 


1.084 


arcsec 


T 


Meridian/transit circle 


26,968 


100.00% 


0.250 


arcsec 


M 


Micrometer 


12,081 


94.56% 


1.742 


arcsec 


X 


Discovery obs. 


9520 


0.81% 


0.898 


arcsec 


H 


Hipparcos obs. 


5494 


100.00% 


0.148 


arcsec 


E 


Occultations 


1736 


100.00% 


0.085 


arcsec 


R 


Ranging 


612 


96.41% 


3.325 


km 


R 


Doppler 


432 


99.07% 


6.660 


km/s 


V 


Roving observer 


372 


49.73% 


0.822 


arcsec 


c 


Encoder 


16 


100.00% 


0.558 


arcsec 


n 


video frame 


12 


100.00% 


0.319 


arcsec 



This improvement can also be illustrated by the evolution of 
the uncertainty on the 2029 B-plane. Table [6] shows the value 
of the nominal position (£, £) and the uncertainty (cr^, o^). The 
uncertainty on the geocentric distance A is also given through 
the value of <x A . One can see that radar measurements enable 
to decrease the uncertainty on the main component £ and on A 
by a factor more than 2, while the uncertainty on £ component 
remains almost unchanged. 



One should note that if the arc data increases thanks to op- 
tical measurements and without radar data, the orbital improve- 
ment using ranging and/or doppler measurements could not be 
that significant. For instance, Apophis has additional optical ob- 
servations from 2011-2012. Considering those data combined 
with the radar measurements, the uncertainty on the distance 
drops to 154 km. If we ignore the radar data, this uncertainty 
value is 177 km. 



Table 6. Uncertainty (o" f ,<r f ) on the 2029 B-plane of Apophis 
considering optical data only or optical and radar data. 







Optical + radar data 


Optical data only 




o> (km) 


6980 ± 15 


6963 ± 17 


(± 


cr ( (km) 


37,440 ± 345 


37,144 ± 788 


A± 


o"a (km) 


38,068 ± 345 


37,791 ± 788 





4.2.2. Future observations 

Gaia is an astrometric mission that will be launched in autumn 
2013. The aim is to have a 3D-map of our Galaxy. There are 
many scientific outcomes from this mission and as far as our 
Solar System is concerned, the satellite will be able to map 
thousands of main belt asteroids (MBAs) and near-Earth ob- 
jects (NEOs) down to magnitude 20. The high precision astrom- 
etry (0.3-5 mas of accuracy on a transit level) will allow orbital 
improvement, mass determination, and a better accuracy in the 
prediction and ephemerides of potentially hazardous asteroids 
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(PHAs). During the 5-year mission, Gaia may observe several 
time the asteroid Apophis (Bancelin et al. 2012b i. Thereafter, 
we performed some tests to quantify the impact of future mea- 
surements (optical, radar or space data) on Apophis orbit when 
the new data are outside the period of observations. Apophis 
will have some favorable conditions of observations in 2013 and 
2021 as it will be at a distance of less than 0.1 au In this context, 
we have considered 5 sets of possible measurements: 

- Si: using all optical and radar data available (2004-2012); 

- Sz- using set Si with additional Gaia data (over the period 
2013-2018) with 5 mas accuracy; 

- S3: using set Si with one additional future radar measure- 
ment in 2013 with l//s accuracy (measurement of a timing 
echo); 

- S4: using set S 1 with one future optical observation done in 
2013 with 0.1 arcsec accuracy; 

- S5: using set S 1 and the case that Gaia would provide only 
one observation with 5 mas accuracy. 

Figure [8] shows that the Gaia data enable to reduce the po- 
sition uncertainty knowledge to the kilometer level (set S2) and 
to keep this value until the close approach of 2029. For com- 
parison, the effect of future accurate measurements (radar and 
optical) can be comparable to the impact of one future Gaia data. 



Set S1 + Gaiadata[5 mas] 
SetS1 + 1 radar data [1 [js] 
Set St + 1 eptidal data [0.1 arcsed] 
Set S1 + 1 Gala data [5 mas] 




Fig. 8. Evolution of the position uncertainty (given by its geo- 
centric distance) of asteroid Apophis considering several differ- 
ent sets of observations. 

In the same context, we have quantified the uncertainty in the 
B-plane of 2029. While the impact of one additional Gaia data 
(set S5) can be compared to the effect of one radar measurement 
(set S3), one set of Gaia observations (set S2) can bring the un- 
certainties around the kilometer level (Tab|7]i. 

Table 7. Uncertainty (<r ( ,cr ( ) on the 2029 B-plane of Apophis 
considering various sets S, of future observations (see text). 





s, 


s 2 


S3 


s 4 


s 5 


<x f (km) 


3.45 


0.23 


3.38 


3.44 


1 


cr ( (km) 


154 


2.7 


10.93 


47 


19.8 



4.3. Impact of Gaia stellar catalogue on orbit uncertainty 

The space mission Gaia will not only provide accurate observa- 
tions of asteroid positions but also and especially a astrometric 
stellar catalogue. The Gaia catalogue (Mig nard et al.|2007) will 
provide unbiased positions of a billion of stars until magnitude 
20 and with an accuracy depending on magnitude (7 fias at <10 



mag; 12-25 yuas at 15 mag and 100-300 /ias at 20 mag). With 
this stellar catalogue, a new astrometry era will be possible. 
New process of reduction of new and archived observations will 
be necessary such as propagation model to third order for star 
proper motions, differential aberration, atmospheric absorption, 
color of stars, etc. With such corrections, an accuracy of 10 
mas is expected on the position of any asteroid deduced from 
Gaia stellar catalogue. In particular, it will be possible to reduce 
future but also old observations with this catalogue even if it 
seems obvious that all the observations for all the asteroids 
could not be reduced again. 

In this section, we studied the impact of this new reduction 
in the orbit uncertainty and we considered the case of (99942) 
Apophis and assumed its current observations: 1518 optical ob- 
servations and 7 radar measurements from 2004 to 2012. 

We assumed that a variable part of the optical observations 
could be reduced again with the Gaia stellar catalogue (i.e. with 
an accuracy of 10 mas systematically) and we computed on fig- 
ure [9j the position uncertainty from 2000 to 2050 depending on 
a percentage of reduced observations with Gaia stellar catalogue 
(0%, 10%, 25%, 50% and 100%). The case of 0% means that 
none of the optical observations could be reduced again and rep- 
resents the current position uncertainty by considering the cur- 
rent accuracy of the observations. 100% means that all the ob- 
servations could be reduced with Gaia stellar catalogue. 



100000 - 




0.01 



2000 



2010 



2020 2030 
date 



2040 



2050 



Fig. 9. Evolution of the position uncertainty of asteroid Apophis 
by considering that a certain percentage of optical observations 
could be reduced with Gaia stellar catalogue, i.e. with an accu- 
racy of 10 mas (see text). 

The variation of the uncertainty is the same for all the 
cases but the amplitude is different (Fig. 15}. With the current 
accuracy of observations, the uncertainty of the Apophis orbit 
is about 10 to 100 km in geocentric distance from 2000 to 
2029 and then sharply increases after the 2029 Earth close 
approach. While the percentage of observations reduced with 
Gaia stellar catalogue increases, the uncertainty of the orbit 
decreases. Figure|9jshows that if only 10% of observations were 
reduced with Gaia stellar catalogue, the accuracy of the position 
would be 20 times smaller than with current observations. If 
the percentage reaches 50%, the position uncertainty would 
be 40 times smaller. Finally, if all the optical observations 
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of Apophis could be reduced again with Gaia catalogue then 
the orbit accuracy would be 50 times better than with current 
observations. 



the quality of the orbit will be greatly improved. For asteroids 
and in particular PHAs, it will help to considerably refine the 
orbit and the impact probability during close approaches. 



In the case where just a few number of observations can be 
reduced again, the choice of the first and last observations ap- 
pears as opportune. Indeed, we have compared this case with a 
random choice of the same number of observations reduced by 
Gaia stellar catalogue. The selection of the first and last obser- 
vations provides the best accuracy. In fact, the first and last ob- 
servations, by virtually increasing the length of orbital arc, bring 
more important constraints on the orbital motion of the asteroid. 

Consequently, we assumed that a small number of the op- 
tical observations could be reduced again with the Gaia stellar 
catalogue and we computed the accuracy of the orbit from 2000 
to 2050. We dealt with five situations: 

- current observations with their current accuracy; 

- we assumed that the first and the last observations can be 
reduced with Gaia stellar catalogue i.e. with a conservative 
accuracy of 10 mas; 

- we assumed that the first five and the last five observations 
can be reduced with Gaia stellar catalogue; 

- we assumed that the first ten and the last ten observations can 
be reduced with Gaia stellar catalogue; 

- we assumed that all observations can be reduced with Gaia 
stellar catalogue. 
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Fig. 10. Evolution of the position uncertainty of asteroid 
Apophis considering that a few number of first and last obser- 
vations could be reduced by Gaia stellar catalogue (see text). 



Figure 10 shows that even if only a few number of first and 



last observations could be reduced by Gaia stellar catalogue, the 
position uncertainty of Apophis can be significantly reduced. In 
particular, if only two observations (first and last) are reduced 
again with the Gaia catalogue, then the orbital accuracy is five 
times better than with the current precision of observations. By 
reducing ten (first five and last five) or twenty (first ten and last 
ten) optical observations of Apophis, the quality of the orbit will 
also be improved. 

The Gaia stellar catalogue will be very useful for astrometry 
because with only few observations reduced with this catalogue, 



5. Conclusion 

Orbital uncertainty parameters can be useful for recovering as- 
teroid, instrument calibration and scientific applications. They 
are also important in order to estimate the conditions of asteroid 
close approaches from Earth and eventualy impact probability. 

In this study, we particularly dealt with the Current 
Ephemeris Uncertainty (CEU) provided by the astorb database 
and U parameter from the mpcorb database. In practice, the CEU 
shows some advantages because it is quickly computable, pre- 
cise and provides a physical value (an angle). Relations between 
CEU, H magnitude, orbital arc and dynamical group have been 
identified. In particular, empirical linear relation between the 
CEU and the orbital arc have been highlighted. The CEU sig- 
nificantly increases while the orbital arc increases when orbital 
arc is between 10 and 250 days and more relatively for 250-8000 
day arcs. For orbital arc larger than 8000 days, the CEU is not 
much improved and reaches about 0.1-0.2 arcsec. For short arc 
asteroids (less than 10 days), the CEU is sometimes hard to com- 
pute and remains large while the orbital arc increases. 

Finally, we present statistical data about astrometric mea- 
surements. Most of asteroid observations (about 94%) come 
from CCD frames. Presently, the more accurate data are radar 
measurements with an accuracy of about 3.3 km in distance and 
6.7 km/s in velocity. Unfortunately, these measurements are still 
rare and mostly limited to PHAs (during close approaches). CCD 
frames remain relatively accurate (0.4 arcsec) and are quite nu- 
merous. These measurements provide the astrometry for most of 
asteroids. The bias in stellar catalogues remains the main cause 
of uncertainty. 

Nowadays, there are two ways for improving orbit uncer- 
tainty. One is to increase the length of orbital arc. We have shown 
how fast CEU decreases while the orbital arc increases. In this 
context, it is necessary to regularly provide asteroid observations 
by performing new observations or dataming. The second way is 
to improve the accuracy of measurements. Radar measurements 
appear as the most accurate but they remain few in number. 

The Gaia mission will be a revolution for asteroid dynamic 
and astrometry. The satellite will provide several tens of obser- 
vations for most of asteroids. For PHAs, in particular the aster- 
oid Apophis, we have shown that Gaia observations will help 
to improve orbit accuracy and to refine its impact probability 
during close approaches. In the B-plane, Gaia will decrease the 
uncertainty at the kilometer level instead of 10 km level cur- 
rently. Moreover, with the Gaia stellar catalogue, it will be pos- 
sible to reduce future and old observations with a very high level 
of accuracy (about 10 mas). Even if only few of the old obser- 
vations could be reduced with this catalogue, the orbit quality 
will be about five times better than currently. In the extreme case 
where all observations can be reduced, the orbit uncertainty will 
be about 50 times smaller than currently. 
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